Pendulum View of The Schrodinger Equation 
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A transformation is found between the one dimensional Schrodinger equation and a pendulum 
problem. It is demonstrated how to construct exact solutions with the resulted pendulum equation. 
The relation of this transformation to the Zakharov-Shabat equations is pointed out. 
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The one dimensional (ID) Schrodinger equation is very 
familiar to physicists, most of whom have learned quan- 
tum mechanics by first solving this equation for some sim- 
ple potentials often used in textbooks 0]. This ID equa- 
tion also finds broad applications in real-world physics 
with some interesting and popular potentials, such as 
harmonic, double- well, and Morse 2] . Additionally, be- 
cause of its underlying mathematical structures, it is 
widely exploited in mathematical physics0,Q. 

I present here a transformation that provides a new 
view for the ID Schrodinger equation. This transfor- 
mation turns the ID equation into a differential equation 
describing a damped and driven pendulum, another well- 
known problem to physicists. Analysis shows that, cor- 
responding to the bound eigenstates of the Schrodinger 
equation, the pendulum equation has a set of special so- 
lutions, which are referred to as critical solutions. The 
number of nodes of a bound eigenstate is related to 
the winding number of its corresponding pendulum mo- 
tion. I also demonstrate how to construct exact solu- 
tions with the pendulum equation, and reveal the rela- 
tion of this transformation to the Zakharov-Shabat (ZS) 
equations Q. 

The one dimensional Schrodinger equation is 



which is called Riccati equation * 



V(x)ip(x) = X 2 if)(x) 



(1) 



where the units are chosen such that 2m = 1 and h = 1 
for simplicity. Writing the eigen-energy E as A 2 is a 
choice of convenience; there is no loss of generality since 
the spectrum of a Schrodinger equation is bounded from 
below. Introduce the transformation, 



ip(x) = \J p(x) sin[a(a;)/2] 



(2) 



where p and a are real. They are different from the am- 
plitude and phase in the Madelung transformation where 
one writes if) = ^fre l{p . To avoid confusion, here p is called 
magnitude and a angle. 

Straightforward calculations show that the transforma- 
tion turns the Schrodinger equation into two equations, 

da(x)/dx = 2 A - 2 A(x) sina(x) , (3) 
dp(x)/dx — 2 A(x) p(x) cosa(x) . (4) 

The function A{x) is related to the potential V(x) as 

V(x) = A 2 (x) - dA(x) jdx , (5) 




FIG. 1: Schematic of a pendulum described by Eq.@. It 
experiences a "gravity" 2A and is driven by a constant force 
2A. 

Three important and interesting observations follow 
immediately. First, Eq.® is self-contained, that is, the 
evolution of the angle a does not depend on the magni- 
tude p. Secondly, with the addition of pd 2 a/dx 2 , Eq.@ 
is in fact an equation describing a driven and damped 
pendulum (see FigQ. For this pendulum, the mass is 
zero (/i = 0), the driving force A is constant, the damp- 
ing comes from a friction proportional to velocity, and the 
fictitious gravity is not constant as depicted by the force 
function A(x). Thirdly, there is a symmetry in Eq.@: 
if a(x) is a solution for a given A, then —a(x) is also a 
solution for —A. As a result, I will concentrate only on 
the case A > from now on. 

For certain types of potentials (e.g., well-shaped), the 
Schrodinger equation can have a discrete set of eigen- 
values (E n = A 2 , n = 0, 1,2, •••). The corresponding 
eigenstates (bound states) are of great importance in 
physics. It is interesting to know what kind of pendulum 
motions these bound states have become as the result of 
transformation (J2J. To find it out, I shall first consider 
one general class of finite force functions A{x) that satisfy 



A{x) < 1 



and 



A(\x\ -> oo) = 1 



(6) 



According to Eq.0, these force functions lead to well- 
shaped potentials that are constant at the infinite bound- 
aries and have "dips" in the intermediate region. More 
general cases, where the function A may diverge at in- 
finities, will be examined later with examples. 
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For this class of force functions, the pendulum expe- 
riences a constant "gravity" at the boundary regimes 
\x\ — > oo, where the equation of motion becomes 

da/dx — 2A — 2 sin a . (7) 

It has fixed point solutions when A < 1. They are 

/ a (A) = 2Nn + arcsin A , /„(A) = {2N + 1)tt - arcsin A , 

(8) 

where N is an integer. All other solutions of Eq.Q ap- 
proach asymptotically to f s . It implies that when A < 1, 
the pendulum motion described by Eq.Q must start 
at a fixed point (x — > — oo) and end at a fixed point 
(x — > oo). There are four possibilities: f s — > / u , / s — > / s , 
/u - * /s, and /„ — > /„. Only the first type are the so- 
lutions (named as critical solutions) that correspond to 
the bound states of the Schrodinger equation since all the 
others have, according to Eq.(@J, divergent p at infinity, 
which violates the boundary conditions. 

The critical pendulum solutions exist only for a dis- 
crete set of A that corresponds to the discrete eigenvalues 
(E n = A 2 )of the Schrodinger equation (JIJ. This point 
can be appreciated by noting that f s is a stable fixed 
point and f u is an unstable one. Given a force function 
A{x) and a driving force A, the pendulum motion start- 
ing at f s almost always comes back, stopping at f s since 
the stable fixed point f s acts as an attractor. In contrast, 
the unstable fixed point /„ behaves like a repulser, where 
the pendulum would be pushed away by any small dis- 
turbance. Therefore, only for some special values A„'s, 
the pendulum is lucky enough to end precisely at f u . In 
this way, the eigenvalue problem of the ID Schrodinger 
equation has become finding the special values A„'s in 
the pendulum equation where the critical solutions, 
fs — ► fu, exist. 

It is well known that the wave function ip(x) is oscil- 
latory, vanishing at positions called nodes. The transfor- 
mation clearly reflects this oscillatory nature. Fur- 
thermore, it establishes a link between the number of 
nodes of the wave function ip(x) and the number of rota- 
tions made in its corresponding pendulum motion: For 
a pendulum making n complete rotations in a given mo- 
tion, the corresponding wave function has n nodes. For- 
mally, the number of rotations is called winding number, 
defined for a pendulum solution a{x) as 

W{\) = {a(x^ oo)-a(a;-> -oo)}/2tt. (9) 

For a given force function A(x), only winding numbers W 
for solutions starting at f s are considered for the sake of 
relevance and interest. As already discussed, for almost 
all values of A, the pendulum in such motions always 
comes back to f s at the end. This means that the wind- 
ing number W(X) is mostly an integer as a function of 
A. Only at the special values A n 's, the winding number 
becomes non-integer. 

Further analysis shows that the winding number W(X) 
should behave as depicted in Fig|2 Starting at zero, the 



winding number is an increasing step function that jumps 
by one at A„'s. Note that the jump at A„ is defined 
as W(A+) — W(X~). This conclusion is drawn on two 
facts. First, as shown in Appendix A, one can prove that 
W(X) an increasing function, that is, W(X') > W(X) 
when A' > A for a given A{x). Since it is integer for 
all values of A except the discrete set, it immediately 
follows that it is an increasing step function. Second, 
we recognize that the jump at A n can only be one as 
dictated by the famous oscillation theorem^, Q. The 
theorem says that the wave function ip n (x) for the nth 
bound states has exactly n nodes. It infers that the node 
number increases only by one between two neighboring 
bound states; consequently, the winding number jumps 
only by one at each step. 
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FIG. 2: Winding numbers of pendulum motions as a function 
of A. Dots denote the winding numbers at A n 's. 

This behavior of the winding number can be used to 
estimate the number of bound states for well-shaped po- 
tentials. Given a force function satisfying Eq.©, one 
integrates numerically or by other means the pendulum 
equation J2J at A = I to find the winding number W, 
which is just the number of bound states for the corre- 
sponding potential V(x). 

One general application of Eq. @ is to construct exact 
solutions. In its phase space shown in Fig|3| the critical 
solutions are represented by curves connecting a pair of 
fixed points, / S (A) and /«(A), for a given A. Accordingly, 
the curves are called critical curves, which are plotted as 
dark lines in Fig|3J The following shows with examples 
how to construct exact solutions from critical curves. 

In one example, a critical curve is chosen for A = 1/ 
Plotted as the solid dark line in Fig|3Ja), the curve is 
mathematically described by da/dx — sin(2 * a — 7r/2). 
After integrating it for a(x), one can calculate out the 
force function A with Eq.@, 

i= ^Q S h(2,)-l (io) 
2 cosh(a;) y cosh(2x) 

According to Eq.©, the related potential is 

_ sinhO) - (V2 + I) sinh(3x) [f - y/2 cosh(2a;)] 2 
2\/2cosh 2 (a;)cosh 3/2 (2x) 4 cosh 2 (a;) cosh(2a;) ' 

(11) 
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FIG. 3: Phase space trajectories of the pendulum as described 
in Eq.© at (a) A = l/\/2; (b) A = 1. Gray lines are for 
constant force functions A while dark lines for a changing 
force function A. Solid circles stand for the stable fixed points, 
squares for the unstable fixed points, and diamonds for the 
fixed points at A = 1. 



which has an eigenvalue at Eq = A 2 , = 1/2. By calcu- 
lating the winding number at A = 1, one knows that 
this is the only discrete eigenvalue. The corresponding 
eigenfunction can be found with Eas. (|3l4|l : it is 



/ [2cosh(2a;)] 1 /2 



(2 cosh x) 



i+a/2 



sm 



(- 



arctan[sinh(2a;)] 



(12) 

In another example, the critical curve is given by 
da/dx = — 2 A + 2 sin a (the dark dashed line in FigEfa)), 
which connects a pair of fixed points, / S (A) and / U (A) for 
any A < 1. The same steps as in the first example lead 
to the potential, 



V(x) = 1 - 2(1 - A 2 )/cosh 2 (Vl - \ 2 x) . (13) 

This potential has only one discrete eigenvalue at A 2 ; it 
is a special case of a well-known potential that has exact 
solutions^]]. Other critical curves, especially, the ones al- 
lowing quadrature integration of da/dx = /(a), can be 
chosen and used to construct potentials for a predeter- 
mined eigenvalue. 

Similarly, one can also construct a potential that has a 
desired number of discrete eigenvalues. The basic strat- 
egy is to use the curves that connect a pair of fixed 
points for A = 1, which are at fi — 2Niir + tt/2 and 



/2 = 27V27T + 7r/2. The potential so constructed will 
have exactly N 2 — N% bound states. The example is 
da/dx = 2\/2sin(a:/2-7r/4) (the dark curve in FigGJb)), 
which connects fi = tt/2 and fa — 5ir/2. With the simi- 
lar steps as in the above two examples, one finds 



V(x) = 



cosh(2\/2a;) - 4 sinh(\/2a;) + 1 
cosh(2V2a;) + 4 v / 2 cosh(\/2a;) + 5 



(14) 



One can verify numerically that this potential has, in- 
deed, only one eigenvalue, approximately at A 2 = 0.5. 

So far, only the cases satisfying Eq.JSJ are considered. 
However, many concepts, such as critical solutions and 
winding numbers, are general and applicable to other 
cases, where A may be divergent at infinities. I demon- 
strate it with an example, A(x) — x, which corresponds 
to a harmonic potential V(x) — x 2 — 1. For this case, the 
pendulum equation becomes 



da/dx = 2 A — 2x sin a . 



(15) 



Consider solutions starting at a = ir as one can check 
that only such solutions have chance to satisfy bound- 
ary conditions for bound states. Through numerical 
calculations, one finds that for almost any value of A, 
these solutions end asymptotically at a = 2(n + 1)tt 
(n = 0, 1, 2, • • •), and the winding number jumps by one 
approximately at A n = y2n (note that E n = 2n is an 
eigenvalue of the harmonic potential V(x) = x 2 — 1). Nu- 
merical calculations show, at these discrete values A n s, 
the pendulum stays close at a — (2n + l)ir for some 
time then strays away, apparently due to numerical er- 
rors. These special pendulum solutions at A„ are exactly 
the critical solutions mentioned before; they correspond 
to the bound states of the Schrodinger equation. De- 
tailed analysis of the asymptotic behavior of Eq. (|15fl in- 
dicates that a — (2n + l)ir acts like a repulser while 
a = 2(n + 1)tt like an attractor even though they are 
not fixed points in general. Finally, look at one special 
case, Ao = 0, where the critical solution is a(x) ~ n, ac- 
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companied by p — e~ x . Combined according to Eq.(|2j), 
they produce exactly the ground state in the harmonic 
potential V(x) — x 2 — 1. 

There is another similar transformation that yields the 
same equations as Eas. (|3l4|) . It is 



ip(x) — p(x) cos[a(x)/2] 



(16) 



However, the potential and the force function are related 
differently 



V(x) = A(x) 2 +dA(x)/dx. 



(17) 



For a given A(x), if the pendulum equation Eq.Q has 
critical solutions at A„'s, then both potentials V(x) and 
V{x) have the same discrete eigenvalues, A 2 's. In short, 
the two potentials V(x) and V(x) have the same eigen- 
values. Note that this conclusion is true only when Eq. 
© has critical solutions for a force function A(x). For 



example, it is not true when A(x) is a step function, for 
which Eq.© has no critical solutions. 

Finally, I point out that the transformations Eas. pilfil) 
are related to the ZS equations In Ref.Q, a trans- 
formation is found reducing the ZS equations to a 
Schrodinger equation; in Ref.yj, another transformation 
is discovered turning the ZS equations into equations sim- 
ilar to Eas. (|3l4f> . Their combination produces these two 
transformations Eas. H2ll6(l . whose generalized versions 
are given in Appendix B. 

The author thanks Junren Shi for helpful discussion. 
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APPENDIX A: WINDING NUMBER W(X) 

For a given force function A(x) , W(X') > W(A) if X' > X. 

Proof: Assume the opposite, W(X') < W(X). Since the 
difference is at least one, the pendulum must make more 
rotations at A than at A'. Note also that at x — > -co one 
has a(x, A') > a(x,X). These two facts infer that there 
must exist a point x, where the two pendulum motion 
functions, a(x, A') and a(x, A) cross each other for the 
first time. At this crossing position, x = x, one should 
have jr{x, A') < 4 s (x, A). However, this contradicts the 
fact, 

dcv nrv 

^(x,X')-^(x,X)=2(X'-X)>0. (Al) 
So, the assumption is not true. We have W(X') > W(X). 



APPENDIX B: TRANSFORMATIONS OF 
ZAKHAROV-SHABAT EQUATIONS 

The Zakharov-Shabat (ZS) equations are 

idUi(x)/dx + (j)(x)U 2 {x) = XUi(x), (Bl) 
idU 2 {x)/dx-(j)*{x)U 1 {x) = -XU 2 (x), (B2) 

where the complex function 4>{x) acts as a potential. We 
choose to write <p(x) = A{x)e lSt ^ x \ where A(x) > and 
S(x) are real. 



Schrodinger equation. - Introduce a transformation 

Ux^e-' 18 ^/ 2 ± i U 2 (x)e tS(x ^ 2 



ip±(x;X) 



^2X + S'{x) 



(B3) 



where S' = dS/dx. By straightforward computation, one 
can verify that ip± satisfy the Schrodinger equation 

- d 2 ^±/dx 2 + V±(x)i/>± = (A + S'/2) 2 i>± , (B4) 

where 

a, AS " 2S^ 3 \X + S')-3(S") 2 



V±(x) = {A Z ±A')± 



2X + S' 



4(2A + 5') 2 



(B5) 

In the special case S(x) = constant, the above equation 
becomes 

- d 2 ip± /dx 2 + V± (x)i>± = A 2 ^± , (B6) 
with much simpler potentials 

V±(x) = A 2 ±A'. (B7) 

These two potentials V± are V and V used in the main 
text. 

Pendulum equation. - For the discrete eigenvalues, the ZS 
equations have \Ui\ = |{7a| = \ff>- This allows one to 
introduce another transformation with 

Ux = Jpe i6 \ U 2 = ^pe ie *. (B8) 

Substituting it into the ZS equations, one obtains 

dp/dx = 2Apsm(9- S) , (B9) 
dO/dx = -2X + 2A cos(0- S), (BIO) 

where 8 = 9i — 9 2 . The above two equations become 
Eas.pi)) when S(x) = with a = tt/2 - 6. 



The combination of Eq. I|B3(I and Eq. (|B8(I results in the 
transformations Eas. i|2ll(j|) . 



[1] L.D. Landau and E.M. Lifshitz, Quantum mechanics 

(Pergamon Press, New York, 1977). 
[2] P.M. Morse, Phys. Rev. 34, 57 (1929). 
[3] Mo-Lin Ge et al, Phys. Rev. A 62, 052110 (2000). 
[4] B.F. Samsonov, Phys. Lett. A263, 274 (1999), B. Bagchi 

and A. Ganguly, math-ph/0302040. 
[5] V.E. Zakharov and A.B. Shabat, Zh. Eksp. Teor. Fiz. 64, 

1627 (1973)[Sov. Phys. JETP 37, 823 (1974)]. 
[6] Einar Hille, Lectures on Ordinary differential Equations 



(Addison- Wesley, Reading, 1969). 
[7] R. Courant and D. Hilbert, Methods of mathematical 

physics (Interscience, New York, 1953). 
[8] W.M. Liu, B. Wu, and Q. Niu, Phys. Rev. Lett. 84, 2294 

(2000). 

[9] B. Wu, J. Liu, and Q. Niu, Phys. Rev. Lett. 88, 034101 
(2002). 



